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Analysis of the Yaw Divergence at P64 Terminus 


INTRODUCTION 


My memo to you on January 12 (Ref„ 1 ) reported on the results of an ex- 
amination of several anomalies and stated that analyses would be made and the 
results published shortly. The analysis of the yaw divergence has been completed 
and is reported here. Analyses of other anomalies about whose causes I was un- 
certain in the preceding memo will be completed and reported separately. 

This analysis is based on rollbacks of a single descent simulation. The 
relative importance of the sources of yaw divergence may be different in other 
simulations. 


SHORT DESCRIPTION OF _G UID AN CE_AND YAWjCONTROL INTER ACTION 

This is intended to provide enough background information on descent 
guidance and control to understand what follows. During lunar descent the guidance 
equations are processed in a ''guidance coordinate frame" (see Fig. 1 ) which ro- 
tates with the moon and whose origin is, on each guidance pass, brought into 
coincidence with the current landing site produced by lunar rotation and landing 
site redesignations, if any. On a nominal descent, the XG, YG, ZG axes of the 
guidance coordinate frame are parallel to the XP, YP, ZP axes of the platform 
frame at the instant of touchdown, but this is true only at that instant because the 
guidance frame rotates with the moon and the platform frame does not. Figure 1, 
adapted from Ref. 2, shows the erection of the guidance coordinate frame. TTF 
is the current time relative to reaching the phase targets (the negative of the time 
to go ) and GAINBRAK = 1 or GAINAPPR = 0 is selected according to phase. Thus 
during the approach phase the guidance coordinate frame is oriented about the 
vertical XG axis such that the YG axis is normal to the vertical plane defined by 
the line of sight to the landing site and the XG axis. The ZG axis is therefore 
horizontal and roughly forward along the direction of motion. 
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Fig. 1 Plan View Showing Orientation of the Guidance Coordinate Frame. 
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The guidance equations produce a window pointing command vector UNWC 
for the powered flight attitude maneuver routine FINDCDUW. For most of the ap- 
proach phase, UNWC is merely the line of sight vector to the landing site. The 
intention is to yaw the LM such that its plane of symmetry — defined by the ZB, 

XB LM body axes - contains the vector UNWC. With UNWC being the line of 
sight vector, the landing site will be superimposed upon the LPD reticle, Line- 
of- sight yaw control works well provided the line of sight is separated from the 
XB axis by a sufficient angle. Figure 2, adapted from Ref. 2, shows why line -of - 
sight yaw control cannot be used all the way to the landing site. As the line of 
sight approaches the XB axis, yaw control becomes indeterminate, and an alternate 
window pointing command vector must be used. The alternate used is the guidance 
coordinate frame ZG axis. The criterion for switching between the line-of-sight 
vector and the ZG axis need not be explained in detail here ( see Ref. 2 ), but for 
a landing which is approximately planar, the line of sight vector is used exclusively 
until the LPD angle reaches 65°, the ZG axis is used exclusively beyond 75°, and 
between 65° and 75° the two vectors are mixed as a linear function of the cosine of 
the LPD angle. 

In a nominal automatic landing, the transition between window pointing 
command vectors starts about 5 seconds before the end of P64, and will just about 
be complete on the final P64 pass. Thus the landing site will be kept on the LPD 
reticle until it disappears out the bottom of the window, and then the LM will yaw 
slightly to aline the reticle in the direction of the ZG axis. The yaw motion pro- 
duced by the final P64 command will normally persist into the second pass of P66. 

Guidance commands to the powered flight attitude maneuver routine 
FINDCDUW consist of a thrust pointing command vector UNFC and the window 
pointing command vector UNWC. Using these two vectors FINDCDUW erects a 
commanded body axis coordinate frame twice, as shown in Fig. 3, adapted from 
Ref. 3. The first iteration satisfies the geometry constraints exactly but fails to 
account for the angular displacement between the thrust vector and the true X body 
axis. The second iteration corrects for this thrust offset and introduces a small 
error in window pointing. This window pointing error, defined as the angle between 
the ZCB XCB plane and the line of sight vector, is the product of the sine of the 
LPD angle and the thrust offset angle about the ZCB axis. Consequently the window 
pointing error ranges from zero when the LPD angle is zero, to the thrust offset 
about Z (whose maximum is under 1 ) w hen the LPD angle is 90°. FINDCDUW does 
not correct this error. 
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Fig. 2 Why Keeping the Landing Site in the Center of Vision cannot be the 
Sole Criterion for Controlling Attitude about the Thrust Axis. 



NOTE: Commanded body axes are computed in two steps. 

Vectors computed the first step are identified by the 
superscript 1. Final values have no superscripts. 


Fig. 3 Geometry of Erection of Commanded Body Axes 
Viewed on a LM Centered Unit Sphere. 
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Using the commanded body axis coordinate frame of Fig. 3, FINDCDUW 
computes corresponding commanded gimbal angles to oring the actual body coor- 
dinate frame into coincidence with commanded frame. FINDCDUW issues to the 
digital autopilot gimbal -angle -increment commands that the autopilot uses to 
increment the desired gimbal angles every tenth of a second during the succeeding 
two seconds. At the end of the two second period the autopilot's desired gimbal 
angles coincide with FINDCDUW' s commanded gimbal angles of the beginning of 
the two second period, and FINDCDUW updates the gimbal -angle -increment com- 
mands from new information. The digital autopilot closes the attitude control 
loop driving the actual gimbal angles into close proximity with its desired gimbal 
angles. 


One statement of the preceding memo (Ref„ 1) was not correct,, The sad- 
den increase in the yaw rate a few seconds prior to P64 was not caused by 
switching from line -of -sight window pointing to ZG-axis window pointing. The 
simultaneity of the break point in the yaw profile with the 65° LPD angle was 
coincidental not causal. The roll angle in the simulation described was caused 
by thrust offset rather than by an out -of -plane thrust pointing command vector. 
Because FINDCDUW does not correct yaw for thrust offset, the yaw attitude was 
not being suppressed by the roll attitude prior to attaining 65° LPD angle as er- 
roneously reported in Ref. 1. 

SOURCES OF YAW DIVERGENCE 


Four sources of yaw divergence have been found: 

1. Out-of-plane velocity due to initial condition dispersions and accelerometer 
bias eventually detected by the landing radar, or due to azimuth landing 
site redesignation. This produces a non-planar approach phase trajectory 
as illustrated in Fig. 1, and the yaw angle acquired is not erroneous but 

is a normal and desirable feature of a non-planar approach. In the run 
analyzed, the Y velocity in guidance coordinates at the start of the approach 
phase was . 31 m/ s to the right and the guidance frame was rotated 4 mr 
about the vertical. 

2. Truncation of the landing site update by the descent guidance equations. 
(This source was discovered by Lowell Hull, Ref. 4. ) Every guidance 
pass the landing site is updated in platform coordinates by the equation 


6 


LAND = LAND UNIT ( LAND + W M X LAND At ) 

where WM is the lunar angular rate vector and At is the guidance period. 
With an approach azimuth north of west, the landing site is updated to the 
right each guidance pass. The updating is truncated to an integer 1/8 m, 
consequently on every guidance pass there is an apparent landing site re- 
designation to the left of up to 1/ 8 m. For a given latitude, the magnitude 
of this truncation error is a sawtooth function of the cosine of the approach 
azimuth (defined to the east of north) as illustrated in Fig. 4. The ap- 
proach azimuth of the Apollo 14 trajectory is 283. 7° and produces a 
truncation error every two seconds of 059 m, about half of the maximum. 
With these repeated apparent landing site redesignations to the left, the 
LM will yaw increasingly to the left, and, regardless of how small the re- 
designation is, the yaw increment will increase each guidance pass and 
become unbounded as the LM flies over the site. Of course, P66 begins 
automatically before this can happen. 

3. The digital autopilot is incapable of attaining zero roll error. Consequently 
any roll error (about the ZB axis) will produce an out-of-plane accelera- 
tion error, rotation of the guidance coordinate frame about the vertical, 
and rotation of the LM in yaw. In simulations, this error has been found 
small compared to the previous two. 

4. A mistake was made twice in the LGC program computations erecting the 
guidance coordinate frame. (See Ref. 5.) The net result of these mis- 
takes is that, with Apollo 14 erasables, on the final pass of P64 the LGC 
uses . 246155 for GAINAPPR instead of 0. This means that the orientation 
of the guidance coordinate frame is based on the out-of-plane velocity on 
the final P64 pass, a violation of the intended procedure. We are favored 
by chance that this gain constant is small compared to unity. 

SIMU LATION RESUL TS 

Simulation results have been analyzed to determine whether or not the 
sources cited explain all of the yaw observed, and we believe they do. Figure 5, 
illustrating these results, contains three parts. The upper part shows that the 
yaw angle and the azimuth angle of the guidance coordinate frame follow very 
closely , as to be expected. The second part of Fig. 5 shows that the landing site 
does move precisely 2. 125 m each 2 seconds, as predicted. The LM motion is 
also plotted, and it converges on the landing site motion, as expected. 


7 



rn oaj d 4 T 7 o* j etji 0 A , a / w f ri/n 

i ‘ : * ' : ; 5 : : 

r ?0 A*\ mnecri i>N 

. p 0 t 6 ~o {ptutie ; c-tcLe 

. 1 ( . 


j 

i ■ : * 

I . ; 

r •••_ 

; 4~ 

L- : - 

-• -j • 'r • 

4- 

- - : - i 

' 1 

• • 

1 

i ; • : i 

h 

i • 



til 

•at 

1 

el 



! . ! 

L : j 

h ; j : 

- . 

- • ; 

; V ... 

; ; 

LA 

van 

_ -...4 


:I,| 

•| 

j. 

1 »■: . 

j’/:? 

... 

; 1 

’ i 

! :• : - 

f - 


.. . 1 

i 







s\ 

A 


3 

£ 

( 

0 


o e,J>tjlx, 

X cfitfXjf. . 
■f dfr* 

o 4 -'5 


Aa/o-lZ 

VE 6 -Ae£J> 


+ ^***-*W+*W 






4 


/t^<f 




© o 


H M X x 


0 oOoo 90o0oo ^ 

xxiK*Xxxxxx*>x 


2,3 

g.a 
2,1 
££ 
/.*? 
/,? 

M 
<6 

iS 

,0>i& 

J>t 40 

,612* 

OOK 

,m<> 

aaJA _ 

,0t™ 

- 

a 


» l 1 1 

I 

- \ • — f 

i- •;• - < -••:•*• •' • • •; 

. • i 

'■ •" ! 1 

; . , 1 

— ♦ — | ' ■ ; ■ \‘"[‘ { 

^^xxxxxxxxxxK^'^^' 

Y M 

/ ■ 

*. UN AA Her to* i*>. a 

; '; i ' ; 

JLF UrC c-*i& Y z sci ^uzs 'jma 

M ">*/•* o « o 0 ! 

PL/WbAM' _ ° •_ _ J 

£00 *}>**/ *' 7 ^ ^ 

; JO •_ ...1 

; i i 

pk Lfrc £*l c Y xn ; Mvh** I* £ ssc 

i • : ■ 1 •■••. : 

' ! : : i 1 

© © © 

i 

' r - T“- ’""• *, ; •■; • ••] ''t ; ; ; j 

( » • ; . . * 

• t ' ....... . , ; 

i ! j ! ; 

i<6' £ O 0 t" '"| 

O ! 

: "' : ' ' ""! ' ' ;■ "' 1 

, ' ■ i i 

p - - 

| j j ; ; * 

. i 1 ’ * • ; ; ; • 

_! j ^ . .;*. .j 4 • { i. 

j i ! j • 1 L ! 


6 - v&A***' f Fn/me 
/{oTMlot^ //^^Afi^eKT 

o* rue <.uaa€t^t r*S£ 


A +i>jj nks 


£ 


T" 


i 

- f- 


* 

qsr> 


«!!?*»»* 

<7?? 7*6 


■ : i rf :X x g o' : 

; ]| 0 0 0 

' if'.?. J 


-It 
. -1 




1 _. 


o coA/fti*** 1 *" 2 rt>\ r_fiyrc*{Ti*** 


OF- C4**9/*6~ Z'T* ufb-*riTF 


a fa otcffi a fr*Tf<rv fvr 

f # f v 4 / i r ? V® - ! »*“ c«w;r<!< l »r>j**o 
A, TW*Cfin>/*. o* LAwJ/f^F i/rf utAAT? 

_ Crji . J . 

/A/ eeaO/ A.*nfi 

A**«AW» cY iC 6g ' 


j LX. ilP 7 VAt, ; AO r%T 7 «A/ W<f-* <t,lwv f 

; i jt/M^MAp^TT L; : TT I' 1 ZT 


t i» | j •• | ‘ I [ I.J. | • 1 1 1 • • j_-^ _j __ 

;+■ JlpW7?«A /A/MfMOl/t <**j AV, 


}. *p.-~ 

i. ... 

177 tPZ. “H® IT* 
Y / A1 ^ -31/000 


if/fH , !TJ 

' t^lwir ; c;< vTlrrt** 






toot. 


/ 0 J+ 


/o*a 



The most interesting revelation of Fig. 5 is the bottom plot which shows 
the rotation increment of the guidance coordinate frame on each guidance cycle 
(X) as compared to a predicted rotation increment (a). The rotation increment 
is within 1 mr of the prediction every pass except the last P64 pass when the 
discrepancy is about 6 mr. The discrepancy on this final pass is due entirely to 
the mistakes in the LGC coding; when the LGC coding is patched to correct the 
effect of the mistakes, the discrepancy becomes zero, i. e. the prediction (D ) 
coincides with the rotation increment ( + ). 

This lower section of Fig. 5 also separates the individual contributions to 
guidance coordinate frame rotation. The circles display the rotation increment 
contributions due to landing site truncation, and the squares represent the com- 
bined effects of landing site truncation and cross -range velocity. In this simula- 
tion, the contributions are about equal except the cross -range velocity becomes 
dominant near the close of P64. The truncation contributions were computed by 
dividing the truncation error by the current range. The cross -range velocity 
contributions were computed by averaging the Y components of velocity in guidance 
coordinates at the start and finish of the current guidance cycle, multiplying by 
the 2 second guidance period, and dividing by the current range. The maximum 
prediction error, (X) - ( a ), is under 1 mr and corresponds to under 1 bit error 
in position. 

The final Figs^ 6A thru 12C, display the yaw response to seven redesigna- 
tion situations for each of three program configurations. The redesignation 
situations are defined on the figures and the three program configurations are as 
follows : 

A. The LGC is patched to cause the apparent GAINAPPR to be close to unity 
on the final P64 pass. This aggravates the program mistake to the maxi- 
mum extent. This was done in such a way as to avoid any other effect on 
the run. For those who are familiar with LGC coding, this effect was 
achieved by replacing JAPFG* +1 by POSMAX after its final use in comput - 
ing TTF/8. 

B. Displays unmodified Apollo 14 behavior. 

C. The LGC erasables are modified to prevent guidance coordinate frame 
erection the final two passes of P64. This is done by loading TCGFBRAK 
with 77776 and TCGFAPPR with ID 14 E+02 B-17 = lO 00257. 
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The dots of Figures 6A thru 12C represent the autopilot's desired yaw 
gimbal angle CDUXD at two second intervals. These figures show that the 
maximum spurious yaw produced by the mistakes in the program using either 
the Apollo 14 erasables or unity apparent GAINAPPR was about 5°. However, 
the behavior under the three configurations was markedly different in every 
case, and the behavior was generally worst with unity apparent GAINAPPR. 


HALVIN G THE MAXIM UM LANDING SITE TRUNCATION ERROR 

It appears that the maximum truncation error could be cut in half by doubl- 
ing the magnitude of the landing site radius ( LAND ) before multiplying by the semi- 
unit vector in the direction of the updated landing site. This would have to be done 
two places; in the computations following TTFINCR and in those following REDES 1. 
In addition it would have to be demonstrated that the redesignation equations could 
never contribute to the truncation error when no redesignation was made, or else, 
if this could not be demonstrated, the redesignation equations could be skipped in 
cases of zero redesignation. 

Considering that the yaw bias seems to work out at about 2° for Apollo 14 
with a . 059 m truncation error per pass, the maximum truncation error of ,125 m 
would probably produce about 4 yaw bias. Is fixing the program worth the effort? 

CONCLUSIONS 


The mistakes found in the program should certainly be corrected for 
Apollo 15 as there is no guarantee we will be as lucky with erasables as we are on 
Apollo 14. Guidance frame erection on the final pass of P64 could be avoided on 
Apollo 14 or 15 by reloading TCGFBRAK with 77776, and there would be no other 
consequence. However, with the Apollo 14 erasables, the consequences of doing 
nothing are benign and that is our recommendation. The maximum landing site 
truncation error could be halved by the minor program change suggested herein, 
but it is doubtful that even this would be worth the effort. All other known sources 
of yaw rotation are normal. 
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Figure 7B 
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Figure 8A 
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Figure 8B 
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Figure 9A 
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Figure 11C 
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